Reentrant transition in the shear viscosity of dilute rigid rod 

dispersions 

Hideki Kobayashj*] and Ryoichi Yamamotclj] 

Department of Chemical Engineering, 
Kyoto University, Kyoto 615-8510, Japan and 
CREST, Japan Science and Technology Agency, Kawaguchi 332-0012, Japan 

(Dated: January 12, 2013) 

Abstract 

The intrinsic viscosity of a dilute dispersion of rigid rods is studied using a recently developed 
direct numerical simulation (DNS) method for particle dispersions. A reentrant transition from 
shear-thinning to the 2nd Newtonian regime is successfully reproduced in the present DNS results 
around a Peclet number Pe = 150, which is in good agreement with our theoretical prediction 
of Pe = 143, at which the dynamical crossover from Brownian to non-Brownian behavior takes 
place in the rotational motion of the rotating rod. The viscosity undershoot is observed in our 
simulations before reaching the 2nd Newtonian regime. The physical mechanisms behind these 
behaviors are analyzed in detail. 
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I. INTRODUCTION 



The viscous properties of dilute dispersions of rigid rods change drastically as the rate 
of applied shear flows 7 increases. Although many previous studies have investigated this 
phenomenon, the mechanism of this viscosity change is not yet completely clear. The aim 
of this paper is to understand the detailed mechanism of the viscosity change by performing 
direct numerical simulations (DNS) for a dilute dispersion of rigid rods that are subject to 
thermal fluctuations in a Newtonian host fluid. 

The relationship between the measurable bulk rheological properties and the microscale 
description of dispersions of rod-like particles has been previously investigated in the litera- 
ture [l-sj]. Giesekus obtained the expression for the bulk stress tensor of diluted spheroidal 
dispersions under shear flow by taking into account the effects of the rotational Brownian 



motion of the spheroids due to thermal fluctuations |3j. Leal and Hinch reported that the 
viscosity behavior is characterized by the aspect ratio I of the rod and dimensionless shear 
rate 'j/D r , where D r is the rotational diffusion constant {5^ 8]. 

In the case of weak shear flow, j/D T <C 1, the dilute rigid rod dispersions exhibit the 1st 
Newtonian behavior, in which the viscosity rj of the dispersion is constant and equal to the 1st 
Newtonian (zero-shear limiting) value 770- For an intermediate regime, 1 <C A f/D T <C / 3 + Z~ 3 , 
the dispersions exhibit shear-thinning behavior, in which 77 oc (7/ Dr)" 1 / 3 . In the case of 
strong shear flow, Z 3 + Z~ 3 <C j/D T , the dispersions reenter the 2nd Newtonian regime, in 
which rj becomes constant again and is equal to the 2nd Newtonian (high-shear limiting) 



value 7700. Similar results have also been obtained in numerical 10|, |ll| and experimental 
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131 ] studies. In this paper, the phrase viscosity transition is used to express the changes 
in viscosity from the 1st Newtonian to the shear-thinning behavior and also from the shear- 
thinning to the 2nd Newtonian behavior. Similar results have been observed for dilute 



dispersions of flexib 



Hinch and Leal 



15| 



5], |f| proposed a theoretical model for the viscosity transitions. They 
considered that the viscosity 77 of the dispersion is determined by the ensemble average of 
the temporal viscosity r)(0,<p) using the probability distribution function (PDF) P^(9,<p) of 
the two orientational angles 9 and </? of the rod, i.e., 

V(i) = J f)(0,tp)P^6,cp)d9dcp. (1) 

Here the form of Pj(9,tp) is shear rate dependent, and the shear rate dependence of the 



dispersion viscosity 77(7) is introduced mainly through this function. 

The rigid rod undergoes a random rotational Brownian motion at low shear rates in the 
1st Newtonian regime, where the effect of thermal fluctuations is dominant over the effect 
of shear flow. Therefore, 

P^(9 , ip) = constant (2) 

holds over the entire phase space of 9 and tp. The viscosity is thus constant with respect to 
the shear rate change in this regime, i.e., 77(7) = 770. 

In contrast, the rigid rod undergoes a deterministic tumbling motion due to strong shear 



flow in the 2nd IN 



fery's equation jl|. Therefore, the PDF approaches the high-shear limiting (non-Brownian) 
asymptotic form with increasing 7, 

p^e,<p) = Pj(d,<p), (3) 

where Pj(9,ip) is the theoretical result jg] derived from Jeffery's equation jl|. The viscosity, 



ewtonian regime. Here the tumbling motion is perfectly described by Jef- 
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therefore, tends to be constant again in this regime, i.e., rj = rj c 

The viscosity exhibits strong shear-thinning behavior in the intermediate regime. The 
PDF is approximately given by 

P^(9, ^ ~ Pj(9, <p) + {D T h)P 1 {9, <p), (4) 

where Pi represents the leading term of the perturbation expansion of the thermal effects. 
It is clearly seen that the contribution from the thermal effects decreases with as the dimen- 
sionless shear rate j/D r increases in this regime, which gives rise to drastic shear-thinning 
behavior. The solid line that is shown in three different flow regimes in Fig. [T] represents a 
schematic illustration of the viscosity transition based on the above considerations. 

Consistent with the theoretical model of Hinch and Leal , the viscosity transition 
from the 1st Newtonian to shear-thinning regimes has already been successfully reproduced 



in various numerical studies [ll|. 



1g| . However, the viscosity transition from the shear- 
thinning to the 2nd Newtonian regime has never been successfully reproduced by numerical 
simulations. For rigid rod dispersions, we could not find any previous studies that have 
been performed at high enough shear rates to approach the 2nd Newtonian regime. Several 
numerical simulations have been conducted for flexible chain dispersions at high shear rates 
that are expected to be in the 2nd Newtonian regime. However, the viscosity transition 



from shear-thinning to the 2nd Newtonian behavior has never been correctly reproduced, 
not even when the hydrodynamic interactions are taken into account using the Rotne-Prager- 
Yamakawa (RPY) tensor 16]. 

In the present study, we used a different class of approach, called the smoothed profile 
method (SPM) 17H21|. that can accurately take into account the thermal fluctuations and 
the hydrodynamic coupling between bead particles with a finite radius a and a Newtonian 
host fluid, based on direct numerical simulations (DNS) of particle dispersions. The viscosity 
of a rigid rod dispersion has been calculated using SPM to reproduce the viscosity transition 
from shear-thinning to the 2nd Newtonian regimes and to examine carefully the validity of 
the theoretical model proposed by Hinch and Leal {5, 6]. 



II. METHODS 



Model 



We solve the dynamics of a single rigid rod in a Newtonian solvent using SPM |17H20|. 
In this method, the boundaries between solid particles and solvents are replaced with a 
continuous interface by assuming a smoothed profile. This simple modification enables us to 
calculate the hydrodynamic interactions both efficiently and accurately, without neglecting 
many-body interactions. The equation governing a solvent with density pf and shear viscosity 
?7f is a modified Navier-Stokes equation: 



Pi 



^ + (u.V)u 



-Vp + r; f V 2 u + fj t <t>{„ + 

fshcar 



(5) 



with the incompressible condition V ■ u = 0, where u(r, t) and p(r, t) are the velocity and 
pressure fields of the solvent, respectively. A smoothed profile function < <f){r,t) < 1 
distinguishes between the fluid and particle domains, yielding <p = 1 in the particle domain 
and 4> = in the fluid domain. These domains are separated by thin interstitial regions, the 
thicknesses of which are given by £. The body force 0f p is introduced to ensure the rigidity of 
the particles and the appropriate non-slip boundary condition at the fluid/particle interface. 



The mathematical expressions for and 0f p are detailed in our previous papers 
The external force f s hear is introduced to maintain a linear shear with a shear rate of 
force is applied with the oblique coordinate transformation based on tensor analysis 
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In the present study, we use a bead-spring model consisting of N freely rotating beads in 
a single rigid rod. The bead diameter is o = 2a. The motion of the ith bead is governed by 
the following Newton-Euler equations of motion with thermal fluctuations: 



M^V^Ff + Ff + Ff + Gr, 



d_ 

7ft 



Rj » i i 



N, H + G- 



(6) 



(7) 



where Rj, V i; and Oj are the position, translational velocity, and rotational velocity of the 
beads, respectively. M; and Ij are the mass and moment of inertia, and F* 1 and are the 
hydrodynamic force and torque exerted by the solvent on the beads, respectively 13, Il8 |. 

and Gp are the random force and torque, respectively, due to thermal fluctuations. The 
temperature of the system is defined such that the long-time diffusive motion of dispersed 
particles reproduces the correct behavior 19|, |20| . Ff represents the potential force due to 



direct inter-bead interactions. 

We use a bead-spring model as a model of rod-like objects with a truncated Lennard- 
Jones potential and a finitely extensible nonlinear elastic (FENE) potential. The truncated 
Lennard- Jones interaction is expressed in terms of U^y. 

12 / \ 6" 



if 





a 



a 



i.i 



e [Tij < 2 e a) 



{nj > 2 e a) 



(8) 



where = |Rj — Rj|. The parameter e characterizes the strength of the interactions, and 
a represents the diameter of the beads. Consecutive beads on a chain are connected by a 
FENE potential of the form 

2" 



^feneO) 



--k c R In 



R 



(9) 



where r = |R«+i — Rj|, k c = 30e/a 2 , and Rq = 1.5a. Ff is the constraint force acting on 
the ith bead due to the bond-angle constraints that cause the connected beads to form a 
straight rod. 

8 N 

FF = 9R-E^-*«)' ( 10 ) 



o=3 

* a = (a - 2)Ri - (a 



DR 2 + R 



(11) 
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where ^ a = is the constraint condition to be satisfied. ii a is a Lagrange multiplier 
associated with the constraints that is chosen such that the condition \& a = is satisfied at 
a time t + h, where h is the time increment of a single simulation step. 

The numerical simulations are performed in three dimensions with periodic boundary 
conditions. The lattice spacing A is taken to be the unit of length. The unit of time is given 
by pfA 2 /?7f, where rjt — 1 and pf = 1. The system size is L x x L y x L z = 32 x 32 x 32. 
The other parameters include the following: a — 4, £ = 2, e = 1, Mj = 47ra 3 /3, iV = 5, 
and h = 6.7 x 10~ 2 . In the presented simulations under shear flow, the Navier-Stokes 
equation is discretized with a de-aliased Fourier spectral scheme in space and with an Euler 
scheme in time 2l| . To follow the motions of the beads, the positions, velocities and angular 
velocities of the beads are integrated with the Adams-Bashforth scheme. The bead particles 
are assumed to be neutrally buoyant, so no gravity effects are considered. At t — 0, the 
rigid rod aligns along the x-axis, which is the flow direction. The total duration r t of each 
simulation is set such that 7T t ~ 3500. The range of k^T is 5.0 x 10~ 4 < k^T < 32 and that 
of 7 is 5.0 x 10~ 3 < 7 < 2.0 x 10~ 2 . From the symmetry of the system, we follow the polar 
angles 9 and tp defined in Fig. [2] to consider the motion of a rigid rod. The angle defined 
between the rod and the x-y plane is denoted by 8, and the angle defined between the rod 
projected on the x-y plane and the x-axis is denoted by ip. 

B. Effective Aspect Ratio 

In the present study, the rigid rod is represented as connected beads. Because the beads 
composing the rod can rotate freely, the effective aspect ratio I differs from the simple 
geometrical aspect ratio L/a, where L ~ Na is the rod's length. Instead, we evaluate / 
numerically with the PDF of the rotating rigid rods without thermal fluctuations in the x-y 
plane, i.e., 9 = 0, as represented by 

= pzi — (12) 

_ sm 

where Co is determined from the normalization condition Jjv Pj(ip)d<p = 1 One obtained 
Eq. (1121) in the following manner. The projection of the PDF Pj(<p) of a rotating rigid rod 
on the x-y plane is governed by a Fokker-Planck equation of the form 

dPM _ djuPM) , n „ apj(y) m , 

6 



where u = (f is the angular velocity of the tumbling rod. When the rigid rod rotates in 
the x-y plane without thermal fluctuations in steady states, the Fokker-Planck equation is 
modified to 

8MW) =0 . (14) 



dtp 



In this case, uj is represented as 



u = y(^TT^ v+ PTi) (15) 

from Jeffery's equation Eq. (fl2|) is obtained because Pj(tp) is inversely proportional to 
u. Figure E] shows that our numerical results of P'((p) = J cos 9P^(9, tp)d9 of the strong 
shear regime agree well with Pj(<p) with / = 7.1. We thus use I = 7.1 for the present rigid 
rod, which is composed of freely rotating beads. 



C. Analytic formula for the viscosity 

Hinch and Leal studied the rheological properties of a dilute dispersion of rigid 

non-spherical particles in steady shear flow. They obtained an analytical formula for the 
dispersion viscosity. We analyze our numerical results with their formula. The dispersion 
viscosity 77 is given by the ensemble average of the temporal viscosity fj(9, (p) using the 
PDF of the two angles for the rotating rigid rods Pj(9, tp), which satisfies the normalization 
condition, 

/ cos6d6 f 2 dipP^(9,ip) = 1. (16) 



The temporal shear viscosity is found to be 

V^,f) = Vf 



1 + $ ( A cos 4 9 sin 2 2<p + 2B cos 2 9 + \- 



H — -F- cos 2 9 sin 2o9 
7 2 ^ 



(17) 



where $ is the volume fraction of suspended particles, D r is the rotational diffusion constant, 

nn 

and A,B,F,I 3 are the shape functions given in previous studies [4H6J]. In the case of rigid 
rod, A, B, F, and ^3 are dependent only on the aspect ratio I. At I — 7.1, A — 8.44, 
B = 0.06, and I 3 = 0.99. 

The shear viscosity of the dispersion is obtained by substituting Eq. ( flT|) into Eq. ([[)). 
When we consider the strong shear case D T <ti 7, we can safely neglect the last term in 



Eq. (ITT]) . The dynamics of the angle ip become decoupled from the angle 9 because the 
angle 9 is sufficiently small for a large I Thus, we obtain the following formula, 



rj(9,ip)P"(9)P'((p)d9d<p 



(18) 



Vi 



1 + $ ( A(cos 4 e) e (sm 2 2(p) v + 2 J B(cos 2 0) fl 
2 



+- 



AE 



(19) 



where (/(0)) fl = J\ cos 9 d9f (9) P" (9), (g(<p)) v = f\d<pg(<p)P'(<p), P"{9) = 

7Y 

Jjv P^(9,(p)d(f, and AE is the error arising from the separation of integrals over 9 and 
tp. We can neglect AE safely because AE is sufficiently small in comparison to the other 
terms. 



III. RESULTS 



In Fig. HI we plotted the intrinsic viscosity 



7/ f $ 



(20) 



of the dispersion obtained from the present simulations as a function of the Peclet number 
Pe. Pe is the dimensionless number that represents the strength of the shear flow normalized 
by that due to thermal fluctuations. In our work, Pe is defined as 

67r?7fcr 3 7 



Pe 



(21) 



We find that the intrinsic viscosity [77] gradually changes from non-Newtonian (shear- 
thinning) to Newtonian behavior with increasing Peclet number, as shown in Fig. HJ The 
present simulation data for [77] show shear-thinning behavior for Pe < 10 2 and 2nd Newto- 
nian behavior for 10 4 < Pe. Those results are in good agreement with previous theoretical 
studies [5, 6, 9]. 

To quantitatively compare our results with those of Hinch and Leal {5, 6], we obtain the 
relation between Pe in our definition and 7/-D r , which is used in Hinch and Leal's work 



instead of Pe. On the basis of the shell model 
for a rigid rod is calculated as 



23 



24l | . the rotational diffusion constant D T 



3(]nl + d(l))k B T 
7rr/fL 3 



(22) 



rf(() = -0.662+ ^1-2^. (23 ) 

In the shell model mentioned above, the contour of the macromolecules of arbitrary shape 
is represented by a shell composed of many identical small beads. The shell model can be 
adequately modeled by decreasing the size of the beads. From Eq. (121"]) and Eq. (122]) . the 
relation between j/D T and Pe is expressed as 

J_ = : p e (24) 

D r 18(ln Z + d(Z)) 1 } 

The theoretical model of Hinch and Leal is plotted also in Fig. |4] with the solid lines in the 
three different regimes, namely, the weak (Ri), intermediate (R2), and strong (R 3 + R 4 ) 
shear regimes. 

According to the work of Hinch and Leal 5| , for the weak shear regime ^ / D T <C 1 , namely, 
Pe <C 7.35 x 10~ 2 , which is denoted by Ri in Fig. HJ [rj\ is constant. On the basis of Ortega's 



work 



251 ]. the intrinsic viscosity [77] of the weak-shear flow regime for a rigid rod with a short 



aspect ratio is calculated as 

'"1 = TE^TTWY (25) 

1.38 8.87 8.82 , . 

T(0 = -0.90 - — + — -. (26) 

This expression is identical to Hinch and Leal's result in the limit of Z — > 00. The lowest 
shear rate that we consider in the present simulations is still not in the weak-shear regime 
because of the extremely long simulation time needed to obtain reliable data. 

For the intermediate shear regime 1 <C j/D z <C Z 3 +Z -3 , namely, 7.35 x 10~ 2 <C Pe <C 26.4, 
which is denoted by R2 in Fig. HI the intrinsic viscosity [77] shows shear-thinning as derived 
from Eq. 

[rj] = dPe" 1 / 3 + 2B + 2// 3 , (27) 

where C\ is an arbitrary constant. Figure H] shows good agreement between the data from 
the present simulation data with that of Hinch and Leal, where C\ was determined to fit 
the simulation data. When I is sufficiently large, the contributions from the last two terms 
in Eq. ( 1271) become negligible, and [77] oc Pe^ 1 ^ 3 . 

For the strong shear regime Z 3 + Z~ 3 i/D T , namely, 26.4 Pe, which is denoted 
by R3 and R4 in Fig. HI the theory predicts that the intrinsic viscosity [77] is constant at 
A(sin 2 2(p} j + 2.B + 2//3 = 3.99 = [77^] from Eq. (fl9l . Here (• • • ) j denotes the ensemble 
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average, which is calculated as 

(/(*>)} J = 1 2 *pf{<p)P 3 {<p). (28) 

We obtained [7700] = 3.82 from our numerical data in the high shear regime, which is denoted 
by R4 in Fig. HJ The error from the theoretical value 3.97 is within 4.26%. 

For the regime 10 2 < Pe < 10 3 , which is denoted by R3 in Fig. [4j the behavior of [7/] 
shows a notable undershoot before reaching the high-shear limiting 2nd Newtonian viscosity. 
This is attributable to the fluctuations in 9. It gives rise to the deviations of P"(9) from 

7T 

its high-shear limiting form P*{9) = f* E Pj(9, (p)d<p and increases with decreasing Pe. The 
rods tend to align in the flow direction with increasing Pe. Therefore, (cos 4 9) and (cos 2 9) 
monotonically increase up to their high-shear limiting values with increasing shear rate. 
This leads to an increase in [77] up to [r]^] through Eq. f fT9|) . To examine the role of thermal 
fluctuations in 9 in more detail, let us define 

M = A(sm 2 cp)j(cos 4 9)e + 2B(cos 2 9)e + T (29) 
= 1.81(cos 4 #) e + 0.12(cos 2 #) e + 2.02 (30) 

to estimate the contribution of 9 fluctuations on the total intrinsic viscosity of the disper- 
sion. Here, P(f) = Pj(f) is assumed in Eq. ( fT9l) . and (cos A 9)g and (cos 2 6% are evaluated 
numerically from the present simulations. The results are plotted in Fig. H] with the square 
symbols. One can see that the data of [rjg] almost perfectly collapse onto those of [77] for 
10 2 < Pe. 

On the other hand, the shear-thinning behavior observed for 10 2 > Pe is attributable 
to the effect of the thermal fluctuations in tp. It gives rise to the deviations of P'((p) from 
its high-shear limiting form Pj(<p) and increases with decreasing Pe. To examine this effect 
quantitatively, we introduce 

[A77] = [77] - M (31) 
= A(cos 4 9) e ((sm 2 2<p) v - (sin 2 2<p)j) (32) 

to eliminate the contribution of 9 fluctuations from the total intrinsic viscosity of the dis- 
persion. Figure [5] shows the behavior of [A 77] as a function of Pe. [At/] decreases with 
increasing Pe, and finally [At/] goes to zero around Pe ~ 150. This value is considerably 
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different from the value Pe = 26.4 predicted by Hinch and Leal [5[ for the viscosity tran- 
sition from the shear-thinning to the 2nd Newtonian but agrees well with our theoretical 
prediction of Pe c = 143, at which the dynamical crossover from Brownian to non-Brownian 
behavior takes place in the rotational motion of the rotating rod at Z = 7.1 2(|. 

IV. DISCUSSION 

Let us discuss the numerical models with which the viscosity transition to the 2nd New- 
tonian regime takes place based on Eq. ffT9]) . For the strong shear regime 1 -C j/D T , we can 
estimate (cos 2 6% ~ 1 and {cos 4 9)e ^ 1 because Pj(9,ip) ~ P*{9)P(tp). Here, P(<p) satisfies 
the Fokker-Planck equation, shown as Eq. ffTBl , for which the formal solution is given by 

P{<p) = d jTdVexp ("4]f /ftM) , (33) 
fty, if) = t/j - (1 - -^j) sin^ cos(^ - 2<p), (34) 

7T 

where C\ is determined from the normalization condition, $2* P{^p)d<p = 1. When A f/D T is 
sufficiently large, P(<p) converges to Pj(<p) represented by Eq. (fT2|) . and the viscosity displays 
2nd Newtonian behavior. 

The above discussion is not valid in the limit of I — > oo, which corresponds to an infinitely- 
long or equivalently infinity-thin rod. In this limit, the angular velocity of the tumbling rod 
becomes zero at tp = from Eq. ( fl5l) . Thus, the rod cannot continue rotational motion 
without thermal fluctuations. This is because the hydrodynamic torque acting on the rod 
becomes zero at if = for I — > oo. Therefore, P(<p) in Eq. ( 133]) is modified to 

Poo(f>) = C 2 J #exp f-^j- [V> - sin-0cos(^ - 2ip)]^J , (35) 

TV 

where C 2 is determined from the normalization condition, Jjv P QO {ip)dip = 1. Using Eq. (1351) . 
the intrinsic viscosity [77] is rewritten as 

M = A (sin 2 2^ + 25 + 2/J 3 , (36) 
(/M)oo= [ 2 dipf^PM. (37) 



It is demonstrated in Fig. [6] that the first term in Eq. ( 136]) shows (sin 2 2</?) tX3 oc (j/D T ) 1//3 for 
the entire range of Pe. Figure [7] shows that only A is increasing with increasing I, while i? and 
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2/1$ tend to be decreasing or constant upon increasing I. We estimate A/(2B + 2/1$) ~ Z L 
for I — > 00. This indicates that 



[r/] oc (sin 2 2^)00 oc (7/A-r 1/3 (38) 

holds for the entire range of Pe without indicating the occurrence of 2nd Newtonian behavior. 
The same conclusion can be derived by considering a characteristic shear rate 7* at which 
the first term in Eq. (136]) becomes comparable to the remaining terms. The condition is 
satisfied at 

£ ~ /5 ' 4 - (39) 

This indicates 7* — > 00 for / — > 00. 



In the case of a previous numerical study 16], the hydrodynamic force acting on each bead 



particle was considered via the RPY tensor. Although the translational hydrodynamic force 
was properly considered, the rotational hydrodynamic torque acting on each bead particle 
was completely ignored in that study. Therefore, it is suspected that the hydrodynamic 
torque acting on the rod becomes zero at if = 0. Therefore, the rod cannot continue 
rotational motion at a high shear rate, where the effect of thermal fluctuations disappears. 
This situation is exactly the same as the case of I — > 00. We expect that the 2nd Newtonian 
regime could be correctly reproduced with the RPY tensor approach if the hydrodynamic 
torque is taken into account properly. 



V. CONCLUSION 



In the present study, we numerically calculated the intrinsic viscosity [r]] of a dilute 
dispersion of rigid rods using a DNS method known as SPM. Simulations were conducted 
under the influence of thermal fluctuations and shear flow in the ranges of 5.0 x 10 -4 < 
k-^T < 32 and 5.0 x 10~ 3 < 7 < 2.0 x 10~ 2 , respectively. We have successfully reproduced the 
viscosity transition from the shear-thinning to the 2nd Newtonian regimes, as was correctly 
predicted by the theoretical model of Hinch and Leal 

There are, however, some discrepancies between the theoretical predictions and the results 
of the present simulations. By defining [Arj\ to eliminate the effects of fluctuations in 9, 
which is not considered in the theoretical model, we confirmed that the viscosity transition 
from the shear-thinning to the 2nd Newtonian takes place around Pe = 150. This value 

12 



is considerably larger than the value of 26.4 predicted by Hinch and Leal |5| but agrees 
well with our theoretical prediction of Pe c = 143, at which the dynamical crossover from 
3rownian to non-Brownian behavior takes place in the rotational motion of the rotating rod 
26]. 

We have analyzed the mechanism of the viscosity undershoot observed in our simulation 
before reaching the 2nd Newtonian regime. Shear flow suppresses fluctuations in ip and 9 
as its rate is increased. The former contributes to decrease [rj\, but the later contributes to 
increase [rj\. The undershoot takes place because of the two competing effects. 

We also conclude that the viscosity transition to the 2nd Newtonian regime can be repro- 
duced correctly only if the hydrodynamic torque is properly taken into account in numerical 
models of the dispersions. 
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FIG. 1: A schematic illustration of the viscosity transition, (a) A typical behavior of intrinsic 

7T 

viscosity [r/] as a function of Pe. Here, P(ip) = J_ 2 vr cos 0P(0, <p)d9 is normalized so that S = 

7T 

J^ K P(ip)d(p = 1. (b) -P(y) in the weak-shear regime where the rod undergoes random tumbling, 
(c) P(tp) in the strong shear regime where the rod undergoes periodic tumbling. 
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FIG. 4: The intrinsic viscosity as a function of Pe. [rj] (circle) and [rjg ] (square). The three solid 
lines correspond to the theoretical result of Hinch and Leal jf], (J: the 1st Newtonian regime denoted 
by Ri, the shear-thinning regime denoted by R2, and the 2nd Newtonian regime denoted by R3 
+ R4. In our simulation, the viscosity shows an undershoot before reaching the 2nd Newtonian 
regime R3. 
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FIG. 5: The behavior of [At]] as a function of Pe. [A77] goes to zero around Pe ~ 150. 
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FIG. 7: The behavior of A, B, and 2/I3 as a function of the aspect ratio I. A (bold solid line), B 
(dashed line), 2/I3 (dotted line). The thin solid line represents I 1 ' 8 . 
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